METHOD OF AND APPARATUS FOR THE ELIMINATION OF THE EFFECTS OF 
INERTIAL INTERFERENCE IN FORCE MEASUREMENT SYSTEMS, INCLUDING 
TOUCH-INPUT COMPUTER AND RELATED DISPLAYS EMPLOYING TOUCH FORCE 
LOCATION MEASUREMENT TECHNIQUES 

The present invention relates to the use of force 
measurement techniques, as in touch screen or other touch 
input computer cathode ray tube (CRT) or related displays 
to locate the points of touch on the face of the display, 
and under adverse circumstances where the display appara- 
tus is rather heavy and imperfectly or poorly or 
insecurely supported or mounted so as to be subject to 
inertial interference effects, such as vibration or sway, 
that may cause errors in the force measurement and result- 
ing computed touch location. 

While directed, accordingly, to obviating such iner- 
tial interference effects and rendering the force measure- 
ments immune to the disturbing effects of movement of the 
mass supported by the force measuring apparatus, as from 
such unsteady mountings or supports for the force measur- 
ing apparatus, the invention even more broadly finds use 
also with analagous problems arising in other force 
measuring applications, as well; such as, for example, 
those of weighing objects on a relatively insecure or 



moving support, or effecting such weighing in the presence 
of a relatively large tare mass or counterweight that is 
excited by rapid operation. ( 

From still a broader viewpoint, the invention relates 
to enabling the accurate measurement of rapidly changing 
forces in situations where the associated mass displace- 
ments can cause significant inertial interference. 

BACKGROUND OF INVENTION 

As above noted, the invention has important applica- 
tion in touch input computer and related displays that 
employ touch force location measurement techniques. In 
particular, it provides significant immunity to inertial 
interference effects in specific systems such as those 
described in my earlier copending U.S. patent application 
Serial No. 857,241, filed January 25, 1992, for METHOD OF 
AND APPARATUS FOR TOUCH - INPUT COMPUTER AND RELATED DIS- 
PLAY EMPLOYING TOUCH FORCE LOCATION EXTERNAL TO THE DIS- 
PLAY. r * 



Before proceeding, however, to address such specific 
applications, it is believed helpful to consider first the 
general problems of the vibration, sway or other movement 
of heavy or poorly or insecurely mounted force measuring 
systems that may introduce errors into the force measure- 
ments and the resulting computations stemming therefrom. 

There is always some amount of mass associated with 
the "measurement side" of a force measuring device to 
which an unknown force is applied, including all or some 
portion of the measuring device itself, and the attachment 
and any related objects between them. If and when the 
device support changes motion, this supported mass is 
generally also accelerated. Since the force creating this 
acceleration passes through the device, the device output 
can no longer accurately reflect the unknown force. 

Additionally, the attachment, the device itself, and 
the device support, deform in some degree when the unknown 
force is applied; usually, in an elastic manner. Thus the 
supported mass must undergo some acceleration as the 
unknown force changes, at which times the device output 
responds not to the unknown force alone, but also to the 
inertial reaction force of the supported mass. 




-4- 

These effects are of practical concern when, for 
example, a force must be measured in the presence of« 
an unsteady support and an appreciable supported mass, or 
a rapidly changing force must be measured in the presence 
of appreciable supported mass and compliance. 

An approach to try to avoid such difficulties may 

SI 

reside in attempting to avoid the combinations of condi- 
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CU tions listed above. Where such combinations are, however, 
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l,d unavoidable, two prior art techniques for amelioration are 

UJ 

r|| noted: viscous damping of the mechanical system itself, 

* ,«% and linear filtering of a measurement signal derived 

^ therefrom. Damping is primarily applicable to situations 

f z [ such as weighing, where a force, once applied, remains 

'its j* 

constant until its measurement is completed. During such 
time as the force remains constant (exclusive of the 
inertial effects)., the intertial error decays exponent- 
ially. When the force must be measured as it changes, on 
the other hand, linear filtering may be particularly 
applicable. The filter may be of analog electronic, digi- 
tal electronic, or even mechanical construction, and is 
typically of low-pass or mixed low-pass and notch design. 
It passes low frequencies unaltered, but blocks those at 



and/or near the resonant frequencies, of the suspended 
mass. A measurement is thus obtained which accurately 
reflects the low frequency part of the unknown force, but 
at the expense of abandoning measurement of the high fre- 
quency part. Both damping and filtering typically provide 
faster measurement than can be acheived with no ameliorat- 
ing technique: but not so fast as could be obtained were 
there no inertial interference present, or if there were 
some way to estimate such interference sufficiently 
accurately so that it could be subtracted from the total 
signal . 

It is to the provision of a novel ameliorating 
approach that the invention is directed, and will be here- 
inafter described in terms of the illustrative system of 
said copending patent application, wherein force measure- 
ment is used to locate finger touches on a typical CRT 
display monitor. Such application necessitates the accu- 
rate measurement of as much as 80 pounds, and must contend 
with flexure in the monitor housing and stand, and the 
table beneath, as well as in the force measuring platform 
itself. It has been found that good results (location 
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errors less than 0.1 inch) can routinely be achieved with 
such a system, provided the table is solid, and the monitor 
housing rests on the platform through small, corner- 
located feet. The lowest resonant frequencies of such 
combinations run in the range of 6 to 10 Hz, and a linear 
filter that strongly blocks frequencies above 5 Hz has 
been found to provide adequate speed of response. 

Such is not, however, the case for less secure moni- 
tor supports. Many monitors, indeed, are supplied with a 
"til t-swivel" stand. This form of support raises the 
monitor center of gravity and softens its support, parti- 
cularly with respect to fore-aft and side-to-side nodding 
motions. These two changes directly affect the two lowest 
frequency eigenmodes, dropping them to about 
3Hz for a large monitor on a tilt-swivel stand. A com- 
parable depression of lowest resonant frequencies is seen 
when a free-standing folding table is used below the plat- 
form and monitor. In addition, disturbances such as bump- 
ing or leaning upon an insecure table can introduce large 
low-frequency inertial errors independent of the free 
resonant frequencies of the system. 



While linear filtering with a lower cutoff frequency 
could, in some instances, be used to eliminate excessive 
interference from these lower frequencies, this would be 
only "at the cost of a touch response delay which is 
unacceptable. On the other hand, with only the higher 
frequency filtering in place, position errors can grow to 
as much as several inches for touches applied when the 
table and/or monitor are shaking from moderate bumping, or 
from several rapidly delivered strong touches. 

As before stated, it should be noted that analogous 
difficulties may arise in other force measuring 
applications, such as that of weighing objects on an 
insecure or moving support, or weighing in the presence of 
a large tare mass which is excited by rapid operation. 

The invention therefore provides a new solution to 
this class of problems. 

OBJECTS OF INVENTION 



It is accordingly an object of the invention, 
therefore, to provide a new and improved method of and 
apparatus for measurement of forces that are immune to 



disturbing effects of an unsteady mounting or support for 
the force-measuring apparatus through canceling or elimin- 
ation of errors caused by such effects in the force 
measurements . 

It is a further object of the invention to provide 
such a novel method of and apparatus for measurement of 
forces that are immune to the disturbing effects of move- 
ment, in response to changes in the applied force, of the 
mass supported by the force-measuring apparatus, as well. 

Another object is to provide these immunities in a 
manner that can be conveniently applied to a wide range of 
circumstances, with characteristics of support and of 
supported mass which are not beforehand known, but for 
which the force-measuring apparatus can be easily field- 
calibrated by the user. 

Still another object is to provide an improved 
touch-input computer or related display employing touch 
force location, including externally to the display, that 
is rendered immune to inertial interference effects. 

Other and further objects will be explained herein- 
after and are more particularly delineated in the appended 
claims . 



SUMMARY OF INVENTION 



In summary, from one of its broader aspects, the 
invention embraces in a touch-input computer and related 
supported display employing touch force location measure- 
ments, a method of eliminating the errors that may be 
introduced into force and/or torque measurements by 
undesired inertial interference motions of one or more of 
the support, mechanical system of the display and/or force 
measuring apparatus itself, that comprises, sensing one or 
more components of force and/or torque applied to the dis- 
play by touch forces to provide force and/or torque 
measurements uncorrected for inertial interference motion 
effects that may arise; sensing lineal and/or rotational 
acceleration of the display in response to such inertial 
interference motions; and correcting the uncorrected force 
and/or torque measurements in response to the acceleration 
sensing to reflect ^substantial elimination from the 
measurements of the effects of such inertial interference. 

Other features of the invention and preferred and 
best mode designs are hereinafter set forth. 
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Underlying the invention, however, are three related 
aspects, corresponding to the above-stated objects. 

1) Particular acce 1 erome t er s are incorporated into 
the force-measuring device in sufficient number and manner 
to record alL relevant degrees of freedom of motion of the 
support (if attached to the support side), or of the 
attachment to the unknown force (if attached to the 
measurement side). A typically different linear combina- 
tion of these ac ce 1 e r ome t e r channels is added in turn to 

! : ! 

!?( each channel of uncorrected force measurement, with co- 

" efficients of combination chosen such that the resulting 
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^ sums reflect the desired force measurements, free of 

error, in the low frequency limit, or in the limit of an 

h= arbitrarily stiff connection to all of the supported mass, 

due to tilting or change of motion of the support. 

2) In addition, or alternatively if desired, sets of 
correction channels are derived corresponding to time 
derivatives of the uncorrected force measurements, includ- 
ing typically at least the set of the second derivatives, 
but possibly also sets of one or more derivatives of other 
orders. A typically different linear combination chosen 



such that the resulting sums reflect the desired force 
measurements, free of errror due to changing the supported 
mass, to the extent that the support itself does not 
change motion. 

3) Coefficients of combination for correction 
channels are determined empirically, in situ. Such 
correction channels may comprise accelerome ter readings, 
and/or time derivatives of various order of these or of 
the uncorrected force readings, and/or related quantities 
of similar value, such as quantities effectively approxi- 
mating the preceding over a desired frequency range, or 
linear combinations of any measured data containing the 
desired acceleration and/or derivative data. The force- 
measuring device is installed into its normal operating 
context, which may vary considerably from installation to 
installation. To calibrate, data are then gathered in the 
same manner as for normal force measurement purposes, but 
while the applied unknown force is allowed to remain at 
zero, and then while the system is intentionally disturbed 
in various ways. The idea is to gather a reasonably rich 
collection of measurement sets from moments when forces 



caused by acceleration of the supported mass are the only 
ones applied to the measurement side of the device, Dis- 
turbances may consist of pushing or shaking the system 
support in differing directions and patterns while data 
are gathered and/or pushing or shaking the supported mass 
at different points and in differing directions; data 
being gathered after the disturbing force is removed. A 
correction matrix may now be found, by methods later 
described, which relates the data collected from the 
correction channels to that gathered from the uncorrected 
force channels. The elements of this matrix comprise the 
desired coefficients of combination to achieve the above- 
stated corrections. 

The invention contemplates that any of these three 
independently novel aspects may be separately used, or 
that they may be used in various combinations. The parti- 
cular preferred application and embodiment described here- 
in benefits" indeed from incorporating all three. In this 
manner, the method and apparatus of my said copending 
application is significantly improved to work accurately 
and reliably with a wide range of massive and softly 



supported computer display monitors. A simple one-time 
calibration procedure is performed at display installation 
time, requiring no equipment, and being effectable in less 
than one minute to complete. For installations that would 
otherwise be plagued by large inertial effects, the typi- 
cal performance improvement achieved is a factor of about 
thirty times reduction in error of reported touch loca- 
tion. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention will now be described with reference to 
the accompanying drawings, 

Fig. 1 of which is a side elevation showing an 

illustrative application in which the preferred 
embodiment supports a computer CRT display monitor 
for purposes, of user touch location as described in 
my said copending patent application; 

Figs. 2a and 2b illustrate the design of an appropriate 

acce ler ome t er for use with the invention, in plan and 
exploded views, respectively; 
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Fig. 3 is an isometric view showing the disposition of 
components within the force-sensing platform, 
employing the ac ce le r ome t e r improvement; 

Fig, 4 is a simplified schematic diagram of electronic 

conversion and calculating circuitry suitable for the 
system of the invention; 

I'll 

'ill Fi S* 5 provides a flow-chart of the firmware operation, of 

m 

;U the preferred embodiment. 
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Fig. 6 illustrates the definition of vector and matrix 
quantities needed to interpert the system in three 
dimensional space, with successive Figs. 6a - e 
showing different components of the measurement 
matrix ; 

Fig. 7 is an equivalent schematic representation model of 
the application system of Fig. 1, including those 
model elements relevant to operation of the preferred 
embodiment; and 

Fig. 8 is a diagram similar to Fig. 7, but with additional 
details allowing for estimation of the accuracy of 
correction . 



supporting table 8 resting upon the floor. Calculations 
performed within 5, transform force and torque measure- 
merits into touch location coordinates, which may be passed 
to a-, controlling computer via cable 9. The controlling 
computer may then respond, via cable 10, with some altera- 
tion of the display image, such as brightening the image 
of a pictured push button whose coordinates have been 
touched . 

Before proceeding to a discussion of the features of 
the invention, a summary orientation of the force-measure- 
ment technique disclosed in my said copending application 
(also incorporated herein by reference) may be useful; 
through the improvement of the present invention is 
separate therefrom and useful also with other types of 
force-measurement systems as before explained. 

Returning, now, to the inertial interference 
suppression, cancellation or elimination concepts of the 
present invention, the stabilizer bar and bracket of Fig. 
1 was provided primarily to prevent overt motion within 
the tilt-swivel base 3, such that once the user has set 
the desired orientation of the monitor, there will be a 
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accurately modeled by a single, dominant, rigid mass, and 
some not; but all high enough in the application of this 
embodiment for linear filtration to suffice for removal of 
their disturbing effects. 

With acceleration readings available and properly 
exploited for motions of bottom plate 7 (or top plate 6), 
it has been found that the corrections relevant to the 
internal motions of the moni t or / p 1 a t f o rm device system are 
mathematically identical to those for a motionless bottom 
plate as later explained. In the abscence of accelero- 
meters, however, motions of the bottom plate 7 and table 8 
must be explicitly considered. Folding tables typically 
allow large side-to-side shearing motion. Some varieties 
allow equally great motion front-to-back. Shelves and 
computer tops often allow large front- to-back nodding 
motion, due to torsion of the support surface. Resonant 
frequencies run 2Hz to AHz and are very easily excited, as 
are non-resonant motions during ongoing disturbance. 

Fig 2a depicts the major elements of a preferred 
accelerome ter of suitable design for the particular pur- 
poses of the present invention, which is shown in exploded 
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assembly in Fig. 2b. A berillium-copper foil 20 is etched 
or punched to the retangular outline shown, with mounting 
holes 21 and slits 22 formed in the same operation. 
Spring action resulting from slits 22 provides region 23 
with mobility perpendicular to the foil plane. Tabs 24, 
thermally isolated by slits 25, provide access for 
soldered or welded electrical connection. A single-sided 
PC board 26 comprises an insulating laminate 27 bearing 
conductive foil region 28, which serves as a capacitive 
counter electrode to region 23 of the free-standing foil 
20. Contact areas 29 are accessible along a thin strip 
not covered by foil 20 after assembly. Elements 20 and 26 
are held in opposition by spacing washers 30 and tubular 
rivets 31, with overtravel protection provided by a 
punched laminate bar 32. Wires 33 provide connection to 
the measuring electronics, and are typically routed some- 
what apart, so as to minimize stray capacitance there be- 
tween. The measuring system chosen, however, is largely 
unaffected by stray capacitance to the grounded enclo- 
sure. The dimensions and thickness of foil 20 are chosen 
to provide a resonance of 23 Hz. This is high enough to 
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provide good fidelity in rendering acceleration below 5 
Hz, yet low enough to avoid excessive aliasing energy 
within 5 Hz of the sampling rate. 

The resonant frequency is associated with a particu- 
lar distance such that region 23 sags under its own weight 
when the assembled accelerorae ter lies horizontally. The 
thickness of washers 30 is chosen to be about half again 
this distance. The area of opposed conductors 23 and 28 
is chosen to give an operating capacitance range directly 
compatible with that of the main sensors (55 opposing 56 
in the upper and lower force-sensing platform castings of 
Fig. 3, later more fully described). The application at 
hand requires acce ler ome te r s with response free of signi- 
ficant amplitude or phase error from about 0.1 Hz to about 
5 Hz. The DC response characteristic of the design chosen 
for this embodiment is very helpful in this regard. 
Although DC coupling and capacitive technique are con- 
venient here, it should be noted that the method of the 
invention may employ acceleration responsive sensors of a 
wide variety of other types and principles of operation, 
if desired. 



-2/- 



Fig. 3 depicts the force-sensing platform 5 seperated 
into its top plate casting 40, and bottom plate casting 
41, each with its associated sub-assemblies. Eight 
springs 42 have pressfit upper and lower still attachments 
43 and 44 respectively. Upper attachments .43 are held to 
bosses 45 by self- threading screws (not shown), with 
screws 46 (one shown) passing through lower casting bosses 
47 into threaded lower attachments 44. When thus 
assembled, and with connector 51 attached to circuit board 
53 through mating connector 52, the platform top plate and 
monitor are supported entirely through springs 42, the 
stiffness of wires 54 and additional supported monitor 
cables, such being insignificant compared to that of 
springs 42. Lower flange 48 surrounds upper flange 49, 
with lip 50 covering what would otherwise be an upward 
opening gap. This arrangement prevents incidental contact 
with keyboard, cables, or other small tabletop objects 
from passing extraneous forces to the top plate 40. Capa- 
citive sensors of displacement (and thereby of force) are 
located in pairs in the corners, each comprising the 
before-identified lower plate 55 which opposes an upper 
plate 56 upon assembly. The capacitor sensor plates 55 



-22- 



Si 
rii 



□ 



and 56 comprise conductive foils on insulating subtrates. 
These sixteen substrates, plus the six accelerometer' 
modules 57, are preferably adhesively mounted to 45-degree 
cast'ing bosses, not shown. 

The six accelerometer modules ^are the important 
elements of Fig. 3 that pertain to the operation of the 
present invention. These are mounted in pairs, as shown, 
with axes of sensitivity inclined 45 degrees to the verti- 
cal, and each 90 degrees to the opposite member of its 
pair. The pairs are positioned to approximate the ver- 
tices of an equilateral triangle, and are turned to make 
the axes of sensitivity lie perpendicular to a line drawn 
from the triangle center. In practice, minor adjustments 
to this pattern may be made to avoid interferences. The 
triangle is made as large as may be accomodated within the 
base, to acheive a. balance relationship between the rota- 
tional and translational sensitivities of the ensemble. 
This arrangement of accelerometer positions and orienta- 
tions is, of course, only one of many, which will be evi- 
dent to one skilled in the art, that share the advantage 
of encoding the six degrees of freedom of rigid accelera- 
tion with excellent linear independence. 
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A simplified schematic of suitable electronic conver- 
sion and calculating means for effectuating the functions 
of the invention is presented in Fig. 4, the circuits 
being carried on circuit board 53 within the platform 
device 5. This is similar to that used in my copending 
application and need not therefore be described in great 
detail. Circuits 85 and the multiplexer inputs 94 have, 
however, been increased in number (to fourteen), to acco- 
modate the six acceleromet er channels 94b and eight dis- 
placement/force channels 94a. Timing circuit 83 has been 
modified to provide somewhat higher frequencies, suchthat 
forty complete conversions of all fourteen channels may be 
completed every second. Finally, microprocessor system 80 
has been enhanced with firmware to develop timeder ivat ive 
correction channels, and to employ accelerometer and time 
derivative correction channels to improve force and touch 
location data made available over communications cable 9, 
Fig. 1. The capacitance 57 in circuit 85, furthermore, 
will be replaced by 55 and 56 for the inputs 94a. 

It is now in order to explain how circuits and compo- 
nents operate to obtain the desired functions. In Fig. 5, 
a flowchart of firmware is presented that may be used in 
the microprocessor system 80 to transform the sensor input 
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data into force and the torque substantially free of iner- 
tial error. The subsequent processing to yield touch 
location coordinates, as described in ray said copending 
application, is not germaine to this operation and there- 
fore,, not repeated in this application. 

The conversion methods of Fig. 4 sample the fourteen 
input channels at fourteen different points in time, 
more or-less evenly spaced throughout the interval of one 
full sample cycle. Within each fourteenth of a cycle, the 
exact moment of measurement depends somewhat on the actual 
value measured. The resulting differential group delay is 
removed by processing step 60. Step 60 may comprise four- 
teen different two-tap finite impulse response ("F.I.R.") 
filters, each performing a linear interpolation of the two 
most recent component values to a common equivalent sample 
time. At the cost of slightly greater total delay, 
three-tap quadratic interpolating filters may be used to 
acheive substantially less differential amplitude 
response . 

In step 61, Fig. 5, by contrast with step 60, the 
same response is employed in each of the fourteen filters. 
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Th i s response is designed to remove frequencies above 
about 5 Hz. In one embodiment, for example, 12 tap 
coefficients form a triangular impulse response, compris- 
ing „six values in descending sequence. This is equivalent 
to successive applications of a six-wide and a seven-wide 
square impulse response filter, which give notches at, for 
example, 5.7 Hz, 6.7 Hz, and integral multiples of these. 
As noted in my said copending application, a linear filter 
which is applied uniformly to all components of the 
applied force and torque vector leaves the computed loca- 
tion of force application unaffected, except for possible 
delay in response. Further, the linear nature of the pro- 
cessing steps in Fig. 5 provide that the uniform applica- 
tion of a given filtering characteristic to all channels 
in step 61 is equivalent to applying the same characteris- 
tic instead to the six channels output from 70, insofar as 
final output 74 is concerned. 

The compensatory differential filtering of step 60 is 
thus provided first to assure that the effective aggregate 
filtering of all channels is identical, and then an add- 
itional common characteristic is introduce in step 61 to 
remove frequencies higher than needed for adequate re- 
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sponse speed, which frequencies may also be too high for 
accurate handling under certain simplifying assumptions 
exploited by later processing steps (such as the use of 
s e c o,q d central difference to approximate the second deri- 
vative.) 

Following step 61, the eight channels responsive to 
Q the force sensors are standardized by linear transforma- 

1 I s 

it\ tion 62, and the six channels responsive to the accelero- 

Lii 

i : f| meters are standardized by linear transformation 63. 

I . t 

j'j'j These transformations are assigned the task of compensat- 

i'ft 

] T ing for mechanical and electronic .component variation from 

one particular platform device to the next, so that if 
multiple devices are used in a set of identically con- 
figured installations, the installation specific calibra- 
tion data may. be identical across the set. Transformation 
62 has an additional task of reducing the eight force sen- 
sor channels to six, in such a manner as to eliminate the 
effect on any of these channels of distortion of the 
mechanical structure of the platform device. The forty 
eight coefficients of combination required by 62, and the 
thirty six required by 63, are determined by a factory 
calibration procedure and stored in a non-volatile memory 
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which is part of the microprocessor system 80. The matrix 
of coefficients for 62 may be described as in my said 
copending application. Determination of the matrix of 
coefficients for 63 is simplified by the fact that the 
only significant source of variation resides in the over- 
all scalar sensitivity of each channel, as opposed to the 
position and direction of its axis of sensitivity. These 
channel sensitivities may then be determined by supporting 
the platform device briefly at each of two fixed angles, 
and noting the differing effect of the acceleration of 
gravity on each channel for the two cases. A diagonal 
matrix comprising the ratios of the expected differences 
to the observed differences may then be used for 63. 
Alternatively, 63 may be the product of this diagonal 
matrix with a fixed matrix describing a transformation to 
familar rectangular coordinates. 

The six components of the most recent previous output 
of 62 are held in delay 64, while the second most recent 
are held in delay 65. Combiner 66 sums, component by com- 
ponent, the current and doubly delayed outputs of 62, 
while subtracting twice the singly delayed output, yield- 
ing second central difference signal 72. Delay 67 is pro- 
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vided to match the group delay of forming the second cen- 
tral difference, such that the signals 72 and 73 remain 
accurately aligned in time. The second central 
difference, before mentioned, is employed as a good low- 
to-medium frequency approximation to the true second deri- 
vative, with the advantage of small group delay. Five or 
seven point approximations may he used in other embodi- 
ments for accuracy to higher frequency, but at the cost of 
greater response delay. The uncorrected force readings, 
available with the appropriate degree of delay from 64, 
then pass to combiner 70, wherein corrections from trans- 
formations 68 and 69 remove the inertial interference sig- 
nal. The coefficients of combination 68 and 69 are deter- 
mined in an installation dependent manner by a user cali- 
bration procedure, performed in conjunction with a soft- 
ware program executed on the controlling computer with 
which the platform device communicates. At the end of 
this procedure, coefficients calculated within the con- 
trolling computer are passed to raicroproccesor system 80 
for storage in its non-volatile memory. 
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The touch locating application for which this embodi- 
ment is intended involves a series of short term measure- 
ments of changing forces, separated by periods of quie- 
scence. An "auto-zeroing" feature 75 can therefore pro- 
vide economy by relaxing the requirement for long term 
stability against measurement offset drift. A variety of 
heuristics may be employed for determining quiescence, 
such as the failure of a number of most recent samples to 
contain any components with a range of variation exceeding 
a certain rainuraum. Quiescent samples identified in this 
way (all of whose components should nominally be zero) are 
accumulated into a set of baseline values. These baseline 
values are presented continuously to combiner 70 for sub- 
traction of static offset from outputs 74, By accumulat- 
ing only a fraction of each quiescent sample into the 
baseline values, a simple recursive filtering effect may 
be obtained which prevents the auto-zeroing logic 75 from 
degrading the signal-to-noise ratio of outputs 74. 



CALIBRATION PROCEDURE 

The calibration procedure for transformations 68 
and 69, Fig. 5, is performed as follows: 

•-. At installation time, the user loads and executes 
a software program supplied with platform device 5 on the 
controlling computer C to which it is attached. This 
program, so operating, is able to pass operating 
instructions and action acknowledgements to the user 
through display device 2. The initial steps of this 
program include establishing communication with platform 
device 5 through cable 9, Fig. 1, and instructing the 
device to enter a calibration mode in which all 
components of signals 71, 72, and 73, Fig. 5, are 
streamed to the controlling computer, rather than the 
usual touch location information. The program then 
monitors the three translational components of signal 71 
for a short while, to establish their quiescent zero 
points. The squared magnitude of applied translational 
force is thereafter obtained from the sum of the squares 
of the deviations of these three components from their 
zero point levels. 

The program next leads the user through a series 
of four to six mechanical disturbances of the system, 
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during which it gathers a matrix of twenty four to thirty 
six measurement vectors, of eighteen components each, 
from which transformations 68 and 69 will be calculated. 
For the first four disturbances, the user 1 is instructed 
to apply finger pressure P at certain approximate 
locations on display monitor 2, Fig. 1. In each case, 
upon detecting an applied translational force of 
appropriate magnitude and stability, the user is prompted 
to remove his finger abruptly, and keep it off until 
instructed to proceed to the next location. During the 
first second after release, six measurement vectors are 
taken from the forty or so available and are retained for 
later processing. The interval between successive 
retained measurements is increased from five to nine 
samples over the interval, to eliminate the chance, for 
some resonance, of strobing only zero crossings. 
Simultaneously, a running sum is kept for the sample 
values of each component during the interval. After the 
six measurement vectors have been gathered, they are each 
diminished by the average vector obtained by dividing the 
running sum vector by the number of samples summed. 

After further fluctuations in the uncorrected 
translational force samples have died away to a low 
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level, the user is instructed to continue with the next 
location. The four locations and directions of finger 
Pressure prior to release are typically as follows: the 
top right corner of the monitor 2 pressing toward the 
rear; the top left corner pressing toward the rear; the 
top right corner, around to the side pressing to the left; 
and the bottom right corner, around to the side pressing 
to the left. However, neither the number nor the 
positions of the locations actually used is critical. 
What is needed is a sufficient number of measurements of 
non-zero free motion of the system, in which each easily 
excited eigenmode (or more precisely, conjugate eigenmode 
pair) makes a number of contributions to different 
measurement vectors at randomly diff ering amplitudes and 
Phases. m addition, the number of data elements 
gathered from correction channels (components of 72 and 
73, Fig. 5 ) should substantially exceed the number of 
coefficients to be determined, so that full constraint of 
the latter may be insured, despite the stochastic nature 
of the independence of eigenmode representation. 

in addition to the twenty four free-motion vectors 
so gathered, an attempt is made to add two more groups 
corresponding to motions driven from the support side 
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The user is instructed to shake the table gently from 
left to right. A group of six measurement vectors is 
collected, as above, except that collection commences as 
soon as the program detects a sufficient level of 
disturbance, and the user is instructed to continue 
shaking through the measurement interval. If sufficient 

t s= - 

|'» disturbance is unavailable with gentle shaking, the 
program will offer no immediate acknowledgement — a 

Til condition to which the user has been instructed to 

Ul 

! 3 i| respond by pressing a key on the controlling computer's 

(I keyboard. Finally, instructions are given to repeat this 

O 

process by shaking from front to back. The driven 
measurements collected in this manner show response to 

("J 

-•J frequencies not constrained to the eigenf requencies of 
free motion; they may be particularly helpful in 
constraining the coefficients of acceleration correction 
matrix 69. 

Were the support soft enough to allow slow, 
substantial motion in all six degrees, the associated 
measurement vectors would together tightly constrain all 
aspects of transformation 69. Few tables of any kind, 
however, allow significant motion in vertical 



-34- 



translation. Similarly, the "visible" (i.e., it affects 
the measurements) excitation of all eigenmodes 
corresponding to the principal supported mass would fully 
constrain transformation 68. Yet the presence of filter 
61 virtually assures that at least some of the higher 
frequency modes will be invisible in the data. 
Nonetheless, effective coefficients for transformations 68 
and 69 can be found from the measurements vectors 
collected as described above, in the sense that they 
correct well for the disturbances actually occurring in 
the installation at hand, while manifestly having no need 
to remove disturbances never found in the data. 

After completion of the measurement-collecting 
Phase, the calibration program enters a calculating 
Phase, during which it finds such coefficients for 68 and 
69. These are then transmitted to platform device 5 for 
storage in the nonvolatile memory of microprocessor 
system 80. 

At this point, the program may proceed to 
determine static calibration coefficients necessary for 
touch force location, as described in my said copending 
application. 
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Consider the equation 

f=f-Af-Ba, (1 
c . u u 

:re." f " and " f 11 denote force measurement six-vectors 
c u 

and 71, corrected and uncorrected for inertial errors, 

n 

spectiijfely; where "a" denotes acceleration six-vector 

?t\ 

; wherie* "A" and "B" denote the six by six constant 

w 

I s : 3 

tricesjj;68 and 69, and where the number of dots placed 
er an ! ihstance of a variable quantity denote the order 
time- 'derivative with which the variable is taken in 

'tar 

e instilnce. Equation 1, then, closely approximates the 

lculat-Eons of Fig. 5 subject to the bandwidth 

Q 

mitatiQn intentionally placed on all variables. 

It has been determined empirically that matrices 

and "B" may be found such that "f " very much more 

c 

osely approximates the- true unknown applied force "f" 

ian does "f ". It is also demonstrated in the 
u 

^reinafter presented analysis, that this is a reasonable 
cpectation. Where, as for the calibration measurement 
.ctors, the applied force "f" and the desired estimate 
" it "f " may both be taken to be zero, we have the 
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Con3ider the equation 

f=f-Af'-Ba, (i 
c u u 

where -f •• and "f " denote force measurement six-vectors 
c u 

74 and 71, corrected and uncorrected for inertial errors, 

respectively; where "a" denotes acceleration six-vector 

73; where "A" and "B" denote the six by six constant 

matrices 68 and 69, and where the number of dots placed 

over an instance of a variable quantity denote the order 

of time derivative with which the variable is taken in 

the instance. Equation l, then, closely approximates the 

calculations of Fig. 5 subject to the bandwidth 

limitation intentionally placed on all variables. 

It has been determined empirically that matrices 

"A" and "B" may be found such that "f very much more 

c 

closely approximates the true unknown applied force "f 

than does "f -. It is also demonstrated in the 
u 

hereinafter presented analysis, that this is a reasonable 

expectation. Where, as for the calibration measurement 

vectors, the applied force "f" and the desired estimate 

of it "f " may both be taken to be zero, we have the 
c 
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but avoids possible problems of overflow due to 
unexpectedly large coefficients. Additionally, should 
modes of disturbance poorly excited during calibration 
ever be strongly excited in practice, the less 
sophisticated determination of »C" may greatly enlarge 
the error, where the alternative at least leaves it 
alone . 



THEORY UNDERLYING THE INVENTION 

To understand how the method of the invention may 
be varied in embodiment, and how it may be applied to 
^ differing force measurement problems, it is necessary to 
examine the underlying theory as part of the teaching. To 
do this, we will examine how lumped parameter mechanical 
models can be used to analyze the touch location 
application, both with regards to the correct functional 
form for inertial correction equations, and with regard to 
the limits on the accuracy improvement which may be 
expected. 

Mechanical schematic diagrams are usefully 
employed, such as Fig. 7^ Fig. 8< r **4-w fr _^. while 
correct equation forms may be readily extracted from 
these by viewing them as one-dimensional physical 
diagrams, they become full six degree representations of 
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the three-dimensional problem with correct interpretation 
of the symbols and associated quantities. Capital 
letters denote six by six real, symmetric, constant 
matrices, while lower case letters denote variable 
six-vectors containing both translational and rotational 
information. The matrices represent not only all 
degrees of those parameters normally associated with the 
component elements (e.g. spring stiffness, mass, moments 
of inertia, magnitude of damping) , but also aspects of 
!•;; relative component position, such that in conjunction with 
these matrices, the diagrams need convey only the 
topology of interconnection to fully specify the 
problem. 

Fig. 6 assists interpretation of vector-valued 
variables and model parameters in three-dimensional 
space. The choice of reference system for these 
quantities is arbitrary, in the sense that the forms of 
the mechanical circuit node equations, and results derived 
from them, are independent of the reference system used. 
Any analysis which requires numerical estimates of the 
component values, however, must choose a convenient 
reference system, and then apply that same one 
consistently to all quantities in the problem. The 
required reference comprises a set of coordinate axes 
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100, Fig. 6a, together with a fixed point R in space. 
Note that axes 100 define not only orthogonal 
translational degrees "1", »2", and "3", but rotational 
degrees "4", "5", and "6", the latter three being the 
outward looking clockwise rotations about the first 
three, in order. This enumeration order is depicted at 
100 as applied to the components of force "f", but 
applies analogously to all six-vector quantities. 

In Fig. 6a, the .set of forces 101 acting on rigid 
body 102 is summarized as to effect, in total state of 
force six-vector »f». An agent at R, supported by rigid 
extension 103 from solid support 104, reports in terms of 
coordinate axes 100, the force and torque components 
required to maintain rigid extension 105 of body 102 in 
equilibrium. The negative of these components, 
enumerated in the order shown at 100, constitute the 
vector "f». (Reversing sign reports the force felt, 
rather than the one applied.) Note that the value of 
"f", while strongly dependent upon the position chosen 
for R, is independent of the path of 105 or its point of 
attachment to 102. Note also that many different force 
patterns 101 may lead to the same value of "f", but that 
the effect of any such force pattern on the dynamics of 
body 102 is completely determined by the value of »f». 
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Fig. 6b shows rigid body 102 displaced from its 
equilibrium p 0sition by displacement six _ vector „, „ 

AH possible rigid extensions of body 102 define a mobile 
but rigid auxill iary frame of reference ^ & ^ ^ 

rigid subsystem 102; or, more briefly, the frame of 102 
At R, the frame of 102 has been marked with a copy 106 of 

coordinate system 100. After ri*«rM a 

Alter displacement, the. first 

three components of "x" maw k~ 

x may be read as the displacement 

of the marked copy of "R" w-t-v, 

^py or r with respect to "r» itself 

Some care i, quired in addressing th£ representation Qf 
rotations; but in the llmlt of small emplacements, the 
ft "' ChrSe COm P°"-ts may be read as that rotation of 
the original coordinates at about each axis in turn 

»hi=h brings the other two axes closest to their marked ' 
copies. Thus the best ali gn in g rotation about the third 
axas develops the component »x ■■, as shown. 

6 

Alternatively, . rotation vector (not shown, may be 
coined which points parallel to the axis of rotation in 
the clockwise direction, with i ength proportlonal to 
angle of rotation about this axis. Components -x », 
"y, and a r e then the projectlons ^ ^ ^ 

this vector on th^ f^cf 

the first, second, and third coordinate 
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axe S/ respectively. Importantly, in the limit of small 
displacements, the sum (component by component) of two' 
displacement six-vectors accurately reflects the 
compounded displacements. 

In Figs. 6c, 6d, and 6e, a mechanical system is 
shown composed of rigid and potentially mass bearing 
subsystems 102, 107, and 108, coupled by massless elastic 
element 109, with characterizing matrix S, and massless 
viscous dissipation element 110, with characterizing 
matrix D. For brevity, elements such as 109 are termed 
springs, and elements such as 110 are termed dashpots; 
these terms imply nothing, however, as to the physical 
form of the elements in question. The spring symbol at 
109, both here and in the schematic diagrams to follow, 
is taken to denote the presence of stiffness in all six 
degrees, not just along the axis of the symbol. Also, it 
may be a summary representation of the effects of a 
number of physically distinct elastic bodies which 
interconnect the rigid subsystems (e.g. 102 and 108) to 
which it attaches. Analogous comments apply to the 
dashpot symbol at 110. Note that while later diagrams 
are purely schematic, Fig. 6 is partly geometric for 
expository reasons. m Figures 6c, 6d, and 6e, spring 
109, dashpot 110, and mass 102 are dissected free of the 
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rest of the system, and connected with massless rigid 
extensions for measurement of their matrix parameters ' 
"S", "D», and "M" as seen at R. 

Returning to Fig. 6c, the agent at R imposes small 
static translations or rotations of 111 with respect to 
112, in each case along one of the six degrees of 
coordinates 100, while at the same time applying whatever 
forces and torques are necessary to preclude motion along 
any of the other five. m each case, the force six-vector 
felt (or change in force six-vector felt, if the initial 
position is not relaxed) , divided by the magnitude of the 
displacement, is entered as that column of the six by six 
matrix »S» corresponding to the number of the component 
displaced. it follows from this definition, then, that 
if "x» is a displacement six-vector of in with respect to 
112, from an initial, equilibrium position, then the force 
six-vector felt through 111 is given by the product 
"S x«. Since »S« is necessarily symmetric, this is also 
the product »x S». A force »-S x» is felt through 112. 

In Fig. 6d, similarly, the agent at R imposes 
three translational velocities and three rates of 
rotation upon 113 with respect to 114. . m Fig. 6e, the 
agent at R imposes three translational accelerations and 
three rotational accelerations upon 115 with respect to 
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116. in both cases, the motions are constrained to one 
coordinate degree at a time, as in Fig. 6c. As before, 
coefficients for the matrix columns of »d» may be 
determined in order, by taking the force six-vector felt 
for each motion (the negative of that applied) , in order 
of degree, and dividing its components by the magnitude 
of the velocity applied. Coefficients for the matrix 
columns »M» may be determined in order, by taking the 
applied force six-vector for each motion, in order of 
degree, and dividing its components by the magnitude of 
the acceleration applied. Like »S»; » D » and »M» are also 
always symmetric matrices. 

"S" and »D» are defined through "felt" forces, 
while »M» is defined through the forces "applied". This 
leads to a sign convention in which the diagonal 
components of - s - and »D« are negative, but also one in 
which terms in node equations of the genera i form 

= M a " 311 conce Ptually of positive sign. Actual 

negative signs appear only in writing terms for which 
the available displacement variable is used in a 
negative sense; that is, it describes the motion of some 
other node with respect to the node of interest. 

in the intact and excited mechanical system, let 
"x" be the six-vector displacement from static 
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equilibrium of 108 with respect to 102, let »y» be the 
six-vector displacement of 102 with respect to 107, and 
let "z" be the six-vector displacement of 102 with 
respect to the inertial frame of R, 104 . Then all of the 
forces acting on rigid subsystem 102 are the force 
six-vector exerted by 109 on 102, which has the value 
"-S x", and the force six-vector exerted by no 

on 102, which has the value »D y». The total force 
acting on 102 necessary to explain the observed 

acceleration of 102 is "M z ". 

Thus the node equation for rigid subsystem 102 
may be written: 

-Sx + Dy = Mz 

* Z (7 

in similar manner, equations for all the nodes of any 
lumped parameter system may be written. These may be 
solved to provide a description of the system dynamics, 
wherein combinations of force six-vectors of interest and 
their time derivatives of various orders are equated to 
combinations of dispacement six-vectors of interest and 
their time derivatives of various orders. 

This completes the analogy between the six-vector 
desciption of six-degree per element mechanical systems 



and the usual scalar valued interpretation of the 
equations of linear system theory, allowing the extensive 
results of the latter to be brought to bear on the more 
compact six-vector notation, (in scalar form, the 
equations are approximately thirty six times as ■ 
voluminous.) 

It is now appropriate to turn to Fig. 7, which 
depicts a lumped parameter model which approximates the 
dynamics of the system of Fig. 1. Platform base 7, the 
top of table 8, and any other objects supported by 8 are 
considered to form a rigid subsystem which is node 134 
(with mass matrix "M^ " ) . Display monitor 2, exclusive of 

its tilt-swivel base, is considered to form a rigid 

subsystem which is node 131 (with matrix "M •■) . The leas 

2 

of table 8 are considered to be an elastic coupling 135 
(with spring matrix "S^") between node 134 and solid 

floor 136. Elastic coupling 133 (with matrix »S ») 

2 

represents the joint effect of the eight springs 42. 

Elastic coupling 132 (with matrix »S ") represents the 

3 

joint effect of tilt-swivel base 3 and stabilizer bar 4 



of Fig. 1. Platform top 6 is considered to form a rigid 
subsystem which is node 139. This is taken for the time 
being as massless, as are bodies forming elastic couplings 
132, 133 and 135. The force »f» to be measured acts upon 
monitor node 131, while interfering force »f acts on 

t 

tabletop subsystem !34. The available data comprise the 
displacement six-vector of subsystem 139 with respect 
to subsystem 134, and optionally, acceleration six-vector 
"a" of subsystem 134. m most actual embodiments, of 
course, available data are invertible memoryless 
transformations of these. 

Massless node 139 may be described by equation 
140, Fig. 7, which leads immediately to expression 141 
for^e^x-vector displacement of node 131 with respect 
to node IsiTyfequttion, ^ 7b an d 7c of nodes !34 and 
131, respectively, may then be written without reference 
to displacement »z». Also, the acceleration ..a" measured 
tor node 134 must correspond to model values as shown in 
Fig-. 7d. 

Manipulation of b*^ 7b and 7c yields ^ 7e h, 
This equates a function of unknown forces only, with a 
Wtion of known displacements only. Combining 7c t 
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and 7d, however, gives ^io. 7 f Thi* ^ 

, . ^ This latter equation has ^ 

.i-Plifx.d dependencies which provide some of the • 
advantages of incorporating accelerometers in accordance 
^th the techniques of the invention. 

Since the uncorrected force readings »f « 

correspond directly to »- S x » >^ 7* 

% x ' 7f may be recast in 

terms of "f » eaa . W^^- 

r as seen m ho. 7o >«,•„ 1 A 
u /g. E\g. 7g corresponds ^1 

<i-t ly t0 the form of previousiy ^^^^ ^ ation 
thus e X p la ini„ g the excellenc r?suits obtained ^ 
Practice in its appl ication to the system of Fig x 
m these equations, the ratio of a second 

-a ..coefficient ratio mat ri,. wlth che general ^ 

• The quantity 3 . describes the ^ 
- • syste m consisti„ g of nothin, but the Mss and spri 

the resonances of the „ass so supported. Where the 
" <*•»<*•>**.<• »y a singie frequency -„., the 
terms « X " and »m s" 1 v'» 

have comparable magnitude when 
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"x" is dominated by components with frequencies near 
"v". For a frequency of one-tenth »v, however, the • 
zeroth order term dominates the second order term by a 
factor, of one hundred (the second derivative scaling as 
the square of frequency, . l„ general, though, the matrix 
has a broad range of eigenfrequencies, and many 
appications, such as touch location, have significant 
sensitivity to error from excitation of ail modes, it 
«y be concluded as a genera! rule, then, that the lower 
order term can be expected to dominate the term of order 
two higher by at least the square of the ratio of the 
lowest eigenfrequency of the coefficient ratio matrix to 
the frequency of interest. As an indication of the 
relative error in one degree of freedom of the final 
system output occasioned by omission of a higher order 
term, this ratio may prove pessimistic by a modest 
factor, since the system is usually designed for good 
sensitivity in the sense that typica! »x" values have 
large projections onto that degree's axis in the vector 
space of »x.., while the error components project 
more-or-less randomly. 

The ratio of uncorrected force terms in >tg. 7g 
(or of displacement terms in H^. 7f, describes the 
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resonances of the mass of 131 Fio 7 ^ 

u±, Fig. 7, with the series 

connection of sp rings 132 and 133 . Interestingly< ^ 
a- ju .t the free mo tions of dispiay 2 in Fig . l ^ ^ 
tabie-is absoiuteiy soiid. The accelerators, as ^ ^ 
the present invention, thus have created , ^ 
Platform ..f^, ,. s ,,, and „ M .. ^ ^^^^ dlsappeared 

fr ° m Vi6W - ^"tio-. t.„ mak es the syste m immU ne 

to inertial interference fro m the support sice, „ hile the 
second order ter m makeS it i™ U ne t0 interference from 
motion induced by changes in the unknown force. 

In a configuration where the lowest of the 
resonances just descried is sufficiency high, the second 
order ter m may b e edited, and the acceleration ter m 

retained to deal with table* , , 

table tilt and motion. Whe n the 

table in reality is solid, of course, the acceleration 
term (and the accelerometers, may be omitted, and the 

second order term retained to deal with a hM 

1 wlcn a heavy monitor 

on a soft base. 

The absence from Via in „ 

om f^. 7 g of constants pertaining 

to the mechanics of the table means that calibration 
-trices exist to correct a given installation for all 
Possible support surfaces, the particularities of the 
installation being l imited to the of ^ 



v 



CO 
CO 

I'M 



ul 
CI 

cn 



display and its relationship to the platform device To 
achieve such a "portable" calibration, hoover, the . 
procedure described earner must be perform on a 
surface mobile enough to allow acceleration in an 
agrees of freedom that might ever plausibly be seen 

Returning to ^. 7e , it can be ^ ^ 

without accelerometers f,,n - 

ecers, full correction can be obtained 

from information i mplc it in the history of the 
uncorrected displacement (equivalents force, data 
"self, subj ect to certain restrictions. The most 
^Portant Imitation resides in the absence of any ability 
to remove error refill * 

or resulting from motion forced through the 

fancies associated with coefficient U3 are quit e 

Also, these equations require 135 to be a true linear 
elastic member, which can be a questionable assumption, 
*»• to stiction with, or bumping against, adjacent walls 
and furniture (the accelerometers mas, this assumption, 

associated with coefficient 1« are low,- thus corrections 
rrom the fourth order derivatives must be made and must 
oe calibrated (which may be done in the same manner as 
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previously described for the oth. 

— also, the caw" ' 

— *us t be modified ^ ; eanS " r system 

— term on the J r 6 S6COnd — 
coefficient 142 "Nation of 

Thls calibration is sn n ^ 
Performed by in situ . StU1 ^"ably 

y m SltU/ empirical means but ma 
additional st S ' but ma y require 

PS Wlth know * non- 2 ero forces wh< k 
require equipment. ' ^ may 

1.1! Further refinements and accural • 

UJ inununity to inertia! ■ * " aChievin * the 

y o inertial interference effects in 

V with the invention may be attained ■ ^ 

§ somewhat more realist , ^ ^ ^ ° f th * 

reaixatxc lumped parameter model of th . 

I dy ° amiCS ° f the ^ual system of Fig 2 rep ^ 
Fi *- 8. Node 13 lof F i g 7h , ^Pr-ntedin 

ut *ig. 7 has been split int-^ 

■»•» -«trl, recr „ ^ 151 ' Node "1 (»ith 

4 "presents the CRT u ^ 

" SU ° h pa »* 0^ the frame and hQus 

treated as c»-„ • n °"s lng as are best 

«~s the POMer supply aM electrMics 
" (along „ ith suoh 

«abil i2er 4 „ are f " me ' h0USi "9' stand 3 and 

^ 1* JtVs treat " " — — 



to describe elastic deformation of the frame and housing 
of the CRT display between 11 and 13. Spring 132 
(Fig. 7) has been replaced by spring 152 (with spring 
matrix-, -g^) in recognition of the fact that spring 132 

of Fig. 7 implicitly incorporated some of the effects of 
spring 153. Accelerometer measurements are presumed to 
be available. 

The node equations, presented ±»-si*. 8 b and 8c, 
lead to force equation E>g. 8 d: ;| ^ V 

Two conclusions are drawn from^rSgT 8d. One is 
that in some applications and embodiments, corrections of 
greater accuracy may be achieved within the method of the 
invention by including corrections from higher order 
derivatives of the uncorrected force measurements, and/or 
from derivatives of the accelerometer measurements, m 
this instance, calculation and calibration may be 
performed exactly as previously described, except for the 
necessary widening of matrix »C» and heightening of »r . 
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The second conclusion is that the errors of the 
Fig. 7 mode! corrected in Fig. 8 are quite minor, and 
that for the actual application of Fig. i, . scheme more 
elaborate than that ot^ 7g is probably not 
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^i" ed - FirSt ' «»«"•* the sub-terms 159 and 160 
<K9- M) that comprise the coefficient of second order 
term 157. Ihe ratio „ f coefficient ^ ^ ^ 

coefficient of term 158 describes the resonances of node 
151 with the series connection of springs 133, 152, and 
153. The same ratio for coefficient sub-term 160 
ascribes the resonances of node 150 with 133 and 152 in 
series. The lowest frequencies associated with 1S9 
then, are much lower than those of l 60; partIy beoa l se 
the scalar mass of 150 is typically at least several 
times smaller than that of l 5 i, and ^ be=ause ^ 
center of gravity of 150 is about three times closer to 
the center of rotation of its effective spring. Lowest 
frequency modes in both cases are primarUy rotations 
about the spring centers, and depend upon moments of 

cavity to spring center distance. The coefficient of 

157 is thus dominated bv 15<3 +-u 

<-gu D y ii> y/ and the ratio Qf the 

-efficient of 156 to this is then close to expression 

'I- Since the mass of 151 is substantially larger than 

t^t of 150 , the ratio of the coefficient of 154 to that 

<* 155 is also approximated by expression 16X. These two 

ratios determine the importance of terms 154 and 156, the 
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two terms that are outside the original form of equation 
1. The mass of the power supply, electronics, etc. i s , 
relatively low, and the stiffness of the housing and 
frame ..is relatively high, especially in the sense of 
there being no axes of spring rotation passing far from 
the center of gravity of 150 to create low frequency 
torsional modes. The lowest frequencies associated with 
expression 161 are therefore found to be in the range of 
tens of Hz. The principal frequencies involved in 
locating even very rapid touches are in the vicinity of 
two to three Hz for the preferred embodiment. Thus 
correction terms 154 and 156 are roughly one hundred fold 
smaller than those already allowed for in equation 1. 

The approach of Fig. 8 may also be used to 
evaluate the effect of bodies other than the power 
supply, to see if the failure of their mass to move 
rigidly with that of "the CRT might require corrections 
outside the original form of equation 1. m the 
application of Fig. !, the top plate 6 is the only other 
candidate for a mass of significant effect. Thus, Fig. 8 
may be reinterpreted with node 150 taken as the mass of 
top plate 6 (Fig. i) instead of power supply 13. 153 is 
taken as base 3 with bar 4, 151 is taken as the bulk of 



display 2, and the equacions „ e sUghciy simpufied ^ 
setting -s" 1 . to 2er0 . For ^ appllcatton ^ ^ ^ , 

«» «» of the top Plate ls about ten times less 

that o f the monitor. (The ratio may be ^ ^ 

a small monitor, bui- a en^i! 

or, out a small monitor suffers less 

interference and requires less 

or the monitor. Toopi-hor *-k 

Together, these imply that the lowest 

CO resonance o f expression l 61 , as here-interpreted, i s 
- ten times higher than that o f the monitor hase 

^nation. 

« - ,t can oe seen that the correction terns tor top 

only about one one-hundredth as much as tha 

ucn as tne correction 
terms present in the original for. of equation ! 

model Di " iPati0n featUr " haVS bee " ° mitt6d '»» ^ 

Z 5yStSmS C ° nSidered " P ° lnt ' «~ W« 0-s 

observed for the lo „ frequency ^ ^ 

greatest correction -i « *-u 

emboH- ^Plication of the preferred 

-odrment are in the range of fifty t0 one hundred 

V"cous dissipation features, however, are readily 
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incorporate, by following the ^ 

- s. This lntroduces terms with ^^^^^^ ^ o ^ • 

order into Che force ^ ^ 

coefficient ratlos to the next iowest ordsr 
express first order ^ ^ 

frequency - the higher tecm may acquire 

;;■ importance to the lower. 

| Progressive^ more elaborate models ^ 

Eli th •*°« r -^ly W9her derivatives, a„ d f0 rce 

« derivatives are li te ly t0 appear on the left side. In 

I CaS "' la " er «» ^nis h ed in a series 

,,, solution formed by repeatedly differentiate 

..p tk '■"■"erentiating and scalina 

the equation, and subtr a( -i-i„„ ^ 

; :=^ ^uuuracting the r^cni+- * 

U „ht«. ■ y_^eresult from itself; 

mSanS " inC °— * — dashpot parallel to 
' ** lnt «"«»« "taxation times (e.g. of the 
Plastic tilt-swivel base taXen by itself J h 

deformed and released, ^ 

released) correspond to frequencies well 
above the range to be corrected Th* , 
of 142 (rt « „ hS l0 " eSt frequencies 

^ 761 ' hOWW «< are much too 1 0 „ for the 
series solution of Sto 7e „ 

?e ( f assumed zero) to have a 
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useful range of convergence. 

in any case, extensions to the original form of 
nation 1 involve the addition of more terms for 
cerivati.es of so,e of t he other orders of the variable 
Six-vectors. These are generally of relatively lo„ 
order, as term significance drops off rapidly with 
increasing order for the f regencies of interest. Series 
-lated terms need usually inc i ude only ^ h ^ 

small to need inclusion. 

When empirical calibration i S used to determine 
6 C ° e " iCiSn " - —ion !, m ori9inal or estended 
J *»». be st fit coefficients win be f0 und given ^ 
;. avarlable, regardless of what model or expectations 

motivated the inclusion of such hc ™, „>. 
^ such terms. Thus 'Eig.. 7 f and 

^- 8d provide somewhat different physical 
interpretations for the coefficient of the second 
derivative of the uncorrected force six-vector 
Progressively more detailed models would provide 
Progressively more elaborate and accurate 
interpretations. Careful spring and mass measurements, 
taxen from progressively finer dissections of the system 
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Of Fig. !, could in prlnciple be used ^ calculate ch . s 
Progression of coefficients. Each such matrix in turn . 
would have a somewhat different numerical value, but 
convening rather rapidly toward the singie ..true., value 
for the second order coefficient. y et this latter is 
just the value found - given that the coefficients of 
other necessary terms are also being evaluated - by the 
-Pirical, i„ sltu cal i b ration method hereinbefore 
described. Determination of what other terms are 
necessary may be approached theoretically, as has been 
seen. zt may also be approached experimentally, by the 
simple expedient of trying forms of equation 1 with 
Offering sets of derivative orders present, to see what 
«r k . best for the class of mechanical system relevant to 
a particular application, it should be noted that 
ommitting an unneeded term of lower order may be as 
important as including a necessary term of higher order 
as ambiguities of solution can otherwise exist which 
Preclude the convenience o, a calibration in which the 
unknown forces applied to the excited system may all be 



zero. 



Finally, m the prefe r re d embodiment of the 
invention, evaluation of the correction terms in 
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equation l i s subiPrf ^ 

sublet to known inaccuracies: a slight 
high frequency roll-up in fh 0 -, 

up xn the accelerometer channels ■ 

reflecting their twenty three Hz r* c « 

y i-iiree hz resonance; the 

variations ln the filtering components ^ ^^^^ 
<> and the use of the second ^^^^ 
- the second derivative. Each of these ^ a source 

^ - - corrections. .Since 

these errors are independent <->«. • 

pendent, their combined expected 

value is around three oerrant 

ee percent, consistent with the before 
-ntioned factor of thirt y improvement attained , 
sediments see*in g improvement of correctlon ^ 

application thus beoin with improving, the accurac y with 
-h t he presMt Mrrection t£rms _ ^ 

than by adding new term. * n 

y new terms to equation 1. 

Therefore, in an alternative, embodiment, samp ling 
I 6 ; 80 " ^ " - —meter resonant 

I:;". " lnCrSaSed ' Md ^ <" ~~ —tors 
each circuit 9l . preferab ly all Similar ratio, ma y 

h .t CreaSed - Ih " e — the accurac. with 

" h 7 h — - °^nal elation 1 are 

-aluated, thereb y further improving overall s y stem 

rejection of inert- -i a i ^ 

j-nertial interference. 
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m stiu « tler embodiment, tlnlng circuit >tfMi ^ 
i- «. circuit of Fl> . , ls (orced to a so _ hat 
amplitude b y factor, test apparatus , and the step 

Ce8? "- SeS ° f C1 " UUS 94 -«or.d. such that cortectl „ s 

lnV " !S reSP °" S " " y »' the filters 60 of 

5. A correcting Inverse over th# ^ ^ 

(up eo abouc cen „„ for acceleronfiter 

"sonence U liaise convoked lnl0 th . accelerometer 

channel filters of 60 Th. n i * 

60. The filters of 60 thus modified 

are convolved into the f liters „, 6,. Ih . second derlva _ 

"7 be. expended .to . „„ e . p flUer 
^e linearity of operations is exploited to change the 
" der " "" — '"'on and matrix multiplication, s0 
"at the second derivative futer , ly be convolued ^ ^ 
«Py of the eight force responsive filters of „. .. „ odi . 
dieted m.tchlng delavs „ ay also be added to the 

OIt " C " annelS - " d «^« Of transformation 60 most 

^ executed to provide both signal 7, and slgnal ?2 q( ■ 

S. Each of the resultent fiiters 6. is then 
*-te„ed in gr o„ p dei ay by truncating the <ver y smaU, 

lead coefficient. I„ this manner corrects 

t»er, correction accuracy is 

improved . 



-63- 



In either of these last embodiments, the low-pass 
cutoff frequency of filters 61 may be raised, „„„„ 
lower group delay. ,„ lhl , nanner _ 

»ey be exploited to provide faster system response. 

»hUe the invention as applied to the .touch-input 
computer and related CRT or similar displays employing 
touch force location measurement techniques has been Illu- 
stratively described, preferably embodying each of the 
novel eccelerometer, derivative correction and calibration 
techniques of the Invention, there are instances where the 
combination of all such techniques may not be required for 
inertia! Interference suppression, as before mentioned. 
Concerning the independence of the accele rome t er and 
derivative connection techniques, there are touch screen 
systems such as some monitor or desktop systems without 
control of the support surface where only the accelero- A 
meter - based correction may be most useful. For a 
general purpose retrofit unit, on the other hand, a major 
Problem resides in large, or softly mounted monitors. Such 
units would benefit from the derivative - based correction 
only, at least for users with a reasonably solid desk or . 
table. The novel calibration procedure is also more 
generally useful. 
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Fr ° m an ° ther «PP««ch or viewpoint, a feature of the 
invention resides in its novel for, of an effective linear 
Altering, at least as applied to force measurement. 
While, as earlier discussed, simple forms of so-called 
linear filtering have been proposed in other types of 
applications, such )U nli ke the present invention, generally 
attenuate all energy near the fluencies of distortion, 
rather than correcting the distortion. The present inven- 
tion, as previously explained, qulte t0 the contraryf 
attenuates the error of the signal but not the signal 
itself, even right at the frequencies of maximum error. 

Conventional filtering, moreover, involves blocking 
frequencies in some range as strongly as possible, while 
Passing others unaffected, and does so over a broad range 
of roughly sim ii ar installations. The method of the pre- 
sent invention, however, requires and involves the genera- 
tion of calibration data that reflects the detailed quan- 
tise characteristics of the supported mass and asso- 
ciated particular mechanical system of the specific 
installation. 

A further striking difference underlying the inven- 
tion is that it requires the output or outputs to have 
generally differing "filter" characteristics with 
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respect to each of the multiple inputs, determined by the 
specific mechanical system characteristics. Only in quit e 
rare instances might two be the same, as previously 
pointed out . 

The novelty of the calibration procedure itself also 
appears to reside in several universal aspects. First, 
there is a large number of numbers to be determined — 
seventy two in the embodiment described — and they must 
be determined all at once, requiring at least the equiva- 
lent of a matrix inverse. A simple matrix inverse, how- 
ever, is generally impossible to use since, in practice, 
there is never enough of the right kind of data to define 
a unique and tailor-made calibration for a specific 
installation. One of many possibilities must thus, under 
the invention, be found which can be expected to perform 
satisfactorily. Finally, the calibration "input" con- 
sists, in accordance with the invention, of approximately 
specified disturbances, with such effected while the real 
system inputs with respect to what is to be measured, are 
allowed to remain at zero. 
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Further modifications will also occur to those 
skilled in the art within the spirit and scope of the 
invention as defined in the appended claims, m add ; tlor 
to application of the several features of the invention 
independently or collectively to avoid the effects of 
inertial interference, these features, as before 
explained, raay be applied to other types of force measure 
ment systems, as well, including, but not lifted to, 
weighing systems and the like. 



